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Abstract 

The applicability of the quasi-static approximation for calculating the two-photon annihilation 
rate of the scalar /q(980) meson envisaged as a KK molecule is critically re-examined. It is shown 
that the validity of this approximation depends on the detailed interplay between the momentum 
dependence of the annihilation amplitude and the momentum space transform of the bound state 
wavefunction of the annihilating pair. The approximation becomes invalid when these two scales of 
variation are similar. An improved method of calculation based on the inclusion of electromagnetic 
corrections to the kernel of the Bethe-Salpeter equation for the interacting KK pair is outlined to 
cover this case and applied to re-evaluate the two-photon decay width for /q(980) in a one boson 
exchange model for the interkaon interaction. The corrections are significant and result in a much 
better agreement with experiment. 
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1. Introduction 



The structure of the lowest mass scalar meson /o(980) has been under debate for some 
time. As possible candidates a qq state [if, a q 2 q 2 state [2], a KK molecule jjjj], or perhaps 
some combination of these structures [ij have been suggested. However, the large two-pion 
decay width predicted for the first option,]? (gg — > tttt) ~ 500 MeV from flux tube-breaking 
(H| to ~ 660 MeV using current algebra jg], does not favor a pure qq configuration in view 
of the typical experimental width [3] of ~ (40 — 100) MeV for the /o(980) — > tth decay. 

On the other hand Barnes jl] , in an earlier paper, advanced cogent arguments in support 
of the molecular picture, using various experimental data to obtain an estimate of the two- 
photon decay width of the /o(980) that agrees qualitatively with experiment when this scalar 
meson is viewed as a KK molecule. As in the case of positronium [9j, Barnes' calculation 
relies on the validity of the quasi-static, or wavefunction at contact approximation, where 
the annihilation rate is computed in Born approximation from the free annihilation cross- 
section multiplied by a kaon current that is proportional to the probability of having a bound 
K + K~ pair at the origin. The latter quantity in turn depends on the assumed two-body 
potential that produces the molecular binding. In [8| this was taken to have a gaussian 
shape, with a depth and range typical of the interkaon Weinstein-Isgur potential [3J . 

Repeating the same calculation for the one boson exchange potential developed in [1] 
one finds that the two-photon width comes out an order of magnitude larger than for 
the Barnes' potential, although the two-pion width calculated in the same quasi-static 
approximation agrees with experiment. We have therefore re-examined the validity of the 
quasi-static approximation by implementing a more complete calculation using the Bethe- 
Salpeter equation. It is shown that the applicability of this approximation depends on the 
detailed interplay between the momentum dependence scales of the annihilation amplitude 
and the momentum space transform of the bound state wave function ip(r) of the annihilating 
pair. If the momentum space transform of ip(r) has a much smaller range of variation 
than that of the annihilation amplitude, the quasi-static approximation becomes applicable. 
In the opposite situation where these two scales of variation are similar, the quasi-static 
approximation can be seriously in error and needs to be replaced by the method of calculation 
outlined in Section 3 below. 



2. The KK bound state 



A non-relativistic approximation for the Bethe-Salpeter (BS) equation for weakly bound 
states T(| was employed in [ij to study the mass and strong decays of the / (980) scalar 



meson considered as a bound kaon-antikaon pair of isospin zero interacting via vector meson 
exchange. In momentum space, this equation reads 



( 



P ' 2 + 2M K - P Q Wp') = * / £E-f (p', p)0(p) (1) 



M K A -uyvvr, 4M 2 J ( 27f) 



where f(p',p) is the two-particle irreducible interaction kernel, or transition amplitude, 
appearing in the BS equation that describes the non-relativistic scattering of the pair in the 
center of mass (CM) system with a momentum change p — > p' for either kaon of mass Mjc- 
The external four-momenta of the BS equation for bound states are off-shell but still 
obey four-momentum conservation. In particular the sum of the initial, equal to final, time 
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components appears explicitly as the energy eigenvalue parameter P on the left hand side of 
Eq. ([1]) and also in the kernel T(p, p') = r(p, p'; Pq)- This eigenvalue lies on the unphysical 
sheet of the complex energy plane where the four-momenta of the colliding particles are off 
their mass shell. It gives the mass and decay half-width P = M — iT/2 = 2Mx + e — iT/2 
of the bound system in the presence of interactions, where e < is the binding energy of 
the KK molecule. In obtaining the non-relativistic form Eq. (pQ) it has to be assumed in 



addition that these four-momenta become "almost" physical for weakly bound states [10J . 

The corresponding eigenvalue equation in coordinate space is given by the Fourier trans- 
form of Eq. (pP). This is just the usual Schrodinger wave equation containing a non-local, 
generally complex, potential with the bound state wavefunction ip{v) having the transform 
4>(p). In [4| an expression for f(p', p) for KK scattering and annihilation has been derived 



using a standard SUv(3) x SUa(3) invariant interaction Lagrangian [TJJ to describe vector 
meson exchange between the kaons in the non-relativistic limit. In coordinate space this 
procedure leads to a real, local one boson exchange potential, plus a pure imaginary contact 
term that describes the KK annihilation into the ttti or nr] isoscalar and isovector channels 
respectively, i.e. to an optical potential of the form 

VopMr) = / dVf (r, r^O « -[l$L(r) + \im(v r a^) 5 3 (r)]^(r). (2) 

Here f (r, r') is the coordinate space representation of f (p',p) and the kaon-antikaon 
annihilation cross section for isospin /. 

Explicit expressions for Vq^ e and in the isoscalar channel that are relevant for the 
present discussion can be found in pj. The resulting optical potential leads to a bound state 
solution for the KK system in an s-state that gives Pq = (981 — 25i) MeV for the mass and 
half width of /o(980) in good agreement with experiment Q. While direct measurements 
of KK scattering and annihilation cross sections have not been reported against which the 
scattering predictions given by V opt can be compared, one can do so indirectly by using two- 
channel unitarity to express the cross section for the inverse process tt + tt~ — > KK, which 
has been measured [12) , in terms of the common inelasticity parameter of the KK and Tin 
channels. This parameter may then be calculated from the S"-matrix given by V op t for the 
isoscalar channel, and leads to a reasonable accord with the experimental cr(7r + 7r~ — > KK) 
cross section 

Taken together, these facts suggest that the optical potential of Eq. fl2]) provides an 
adequate representation of the low energy dynamics of the KK system, at least for / = 0. 



3. Two photon decay of / (980) 

We now revisit the two-photon decay width problem for / (980). One can first implement 
the quasi-static approximation by replacing the v r a^> isoscalar cross section for annhilation 
into two pions with that for K + K~ annihilation into two photons and then calculating the 
resulting imaginary part — iT{ , yy)/2 of the energy shift given by Eq. (j2J) in perturbation 
theory. The result is 0, [1] 

r(77)«r, s ( 77 ) = ^ +jr -(o), v*+a--(o) = 4=^(0). (3) 
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with a = e 2 /47r, after projecting out the charged channel amplitude given by the overlap 
C =< K + K-\KK >= l/y/2 for good isospin from the KK isoscalar ground state wave 
function ip(0) at the origin. In obtaining this estimate we have ignored a possible K°K° — » 27 
contribution that lies beyond the scope of the present calculations that assume point kaons, 
but this contribution is known to be negligible at low energies in any event 0. Using 
Eq. ([3]) one finds values for the annihilation width of ~ 0.6 to 0.9 keV for the Barnes' 
potential depending on how ip(0) is estimated, or the much larger value 5.59 keV for the 
one-boson exchange potential. The experimental value is listed as 0.31 l^ii keV [7|. 

This spread of calculated values simply reflects the impact of the details of the assumed 
KK interaction potential on the value of ip(0) used in Eq. ([3]). The more important question 
is how well Eq. ([3]) actually meets the quasi-static conditions used in its derivation, i.e. that 
of being able to calculate the annihilation while ignoring binding and vice versa, at the 
typical binding energies [7;] of about 10 to 20 MeV encountered for the KK molecule. In 
order to address this question we include the contribution of electromagnetic interactions 
between the charged kaons present in the / (980) ground state pi| as part of the interaction 
kernel of the BS equation by replacing r(p', p) with r(p', p) + r em (p', p) in Eq. ([T]). Since 
the /o(980) — > 27 decay width is related to ImT em , we require two-particle diagrams with 
at least two intermediate state photons. For K + K~ scattering there are five diagrams of 
this type to order a 2 that together form a gauge invariant set. They are shown in Fig. 1. 

Call iTem the sum of these five diagrams. Then one can show after some calculation 
using the Cutkosky rules [Io| that the imaginary part of the electromagnetic contribution 
to Eq. ([1]) may be expressed in a form reminiscent of an optical theorem as 

2,:/ ™/|^ r "" (p ' pmp) = ^/^ M " (p,) /|^ M - (P)W(P) 

= 3^/>E^,(p)/|^M-(p)«(p) (4) 

where M M<T (p) is the tensor amplitude for the annihilation of K + K~ into two photons at 
kaon momentum p; C = 1/V2 is the isospin projection factor as before. The symmetry 
factor 1/2 arises naturally to restrict the integration du over the scattering direction of one 
of the photons to half the solid angle. The second form involving the sum over the photon 
polarization vectors e^(k), follows upon introducing their completeness relation and setting 

M AA ,(p) = [(£(k)M^(pK,(-k)] (5) 

We now obtain the electromagnetic decay width f — > 27 by including Jmf em (p',p) 
in the right hand side of Eq. ([T]) for the charged channel and calculating the additional 
imaginary shift in the total energy P — ► P — iT('yy)/2 that this produces in perturbation 
theory. Then 

r <^> = l§fWf inmimtm[v '' p)Mp) 

^6sW5l/(0« p)M "' (p)|2; (6) 

Note in passing that this form reduces to the quasi-static result in Eq. (j3J) if only the static 
part p — > of the annihilation amplitude is kept. For then M\y factors out of the momentum 
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integral to reproduce the free annihilation cross section as v T a em = J2xx> \M\x\ 2 / (dAirM^) , 
while the remaining integral squared over ((f>(p) just equals ip 2 K+K -{Qi). 

Since the Fourier transforms 0(p) = 4>{p) are spherically symmetric in momentum space 
for s-waves it is convenient to split the integral in Eq. §6§ into radial and angular parts, 

d^p 1 

,^0(p)^AA'(p) = T^-.g , 



1 XX 

where 



©av(p) = JdQM xy (p,Q) (8) 

The transition amplitudes M\\>(p) in these expressions are required at off-shell values 
of the kaon four-momenta. We take these as [^Po, ±p] ~ [M K , ±p] in the CM system after 
ignoring the small binding energy of the kaon pair. The outgoing photons are of course still 
on-shell with four-momenta [Mr-, ±k] where |k| = Mx to the same approximation; each 
photon carries away half the total mass of the decaying bound state. We evaluate M AA /(p) 
at these off-shell four-momentum values from the standard expression for the annihilation 
amplitude M^ a (p) for charged point-like bosons, without form factors, that can be found in 
[13] for example. Then from Eq. ([5]) 

MMP) = e 2 {2(e A . .„) + ^f^^ + < k - "<<>} <») 

where |k| = Mr, after introducing the transverse gauge e^(k) = [0, e\] with k ■ = and 
e A ■ e A' = <^aa' for the photon polarization vectors. 

The angular integral in Eq. (jSj) over all kaon momentum directions is easily performed 
relative to a set of axes where the direction of the photon momentum k defines the polar 
axis and the two remaining orthogonal directions by the two polarization vectors of either 
photon. Then 0aa' (p) has a common diagonal value in the polarization indices and is given 
by0 



Qxxip) = 8ne\,(p)6 xy , 7 >) = [t±£- In - \p 2 \ (10) 



after expressing the three-momentum p in units of Mr-. Thus the final integral I\y that is 
required for the width calculation is also diagonal, with a common value I that we factor as 

roc 

I = 8ira*l){0)R, R= dpf(p)d 11 (p), (11) 

" 

with f(p) defined as the normalized Fourier transform 

Tp 1 /"Co 

'«=4** /. * /(p)=1 - (12) 

As before ip(0) is the radial wave function of the isoscalar KK pair at the origin. Thus 
R = 1 if 9 71 (p) is replaced by its static limit value # 77 (0) = 1. The final form of Eq. (J6j) for 
the two-photon decay width then reads 



2vra 2 



r(rr) = ^n+K-QW = r 9S ( 77 )^ (is) 



The coefficient R 2 directly measures the deviation of this width from the quasi-static limit 
given by Eq. (j3J). 
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4. Calculations for two different potentials 

Let us now recalculate the two-photon decay width for both the one-boson exchange 
(OBE) potential and Barnes' gaussian potential. For computational convenience it is useful 
to replace both potentials by their Bargmann equivalent potentials [13] that have the same 
scattering length a and effective range r as calculated numerically for the original potential 
forms. Then (a , r ) translate into the two parameters (a, b) of the equivalent Bargmann 
potential form given in the Appendix. Analytic expressions for both the bound state wave 
function as well as its Fourier transform are then available in closed form, see Eqs. (fT9|) . (I2T]) 
and Fig. [31 

(i) One Boson Exchange Potential. For this case the scattering length and effective range 
were calculated in [4| . They are ao = 5.835M^ and r = 1.187M^\ giving a = — 0.1936M^- 
and b = 1.491M^. For this parameter set there is a single bound state at —18.60 MeV. 
Then Eq. gives ^(0) = 0.260M^ /2 , that leads to T qs {^) = 5.59 keV. Taking <f>(p) from 
Eq. (121 p to do the integral for R numerically, one finds R = 0.287 so that 

r( 77 ) = 0.46l°i° keV, (14) 

or a reduction in the quasi-static value by more than an order of magnitude. The error bars 
have been estimated by considering the experimental uncertainties [7J on Pq = [(980 ± 10) — 
« (20 to 50)] MeV for the /o(980) mass and 2n decay half width that change the values of a 
and b accordingly 

(ii) Gaussian Potential. In [8| Barnes estimated ^(O) using a gaussian potential form 
V(r) = — V exp(— r 2 /2r|) with V = 440 MeV fitted to the typical interkaon potential 
depth of the Weinstein-Isgur model, and a range of r g = 0.57 fm = 1.435Mj^ in order 
to reproduce a binding energy of ~ —10 MeV for the KK pair. The scattering length 
and effective range for this potential are nearly double those of the OBE potential: one 
calculates that a = 8.404M^ 1 , r = 2.397 Mj^ 1 giving a = -O.U37M K , b = 0.6910M^ 
for the parameters of its Bargmann equivalent. The larger scattering length translates 
into a weaker binding energy, —10.2 MeV, and thus a smaller wave function at the origin, 
V> p (0) = 0.102il4 /2 < ^gV(O). This gives the much smaller value 0.86 keV for T gs (77). 
Doing the R integral as under (i) the result is R = 0.632 so that now 

T( 77 ) = 0.34 keV (15) 

This is only about 2| times smaller than r gs (77). The quasi-static limit is thus more 
reliable for kaons moving in the Barnes' gaussian potential than in the one boson exchange 
potential. The reason for this is clear from Fig. [2j This figure compares the normalized 
Fourier transforms f(p) of Eq. (fl2"l) for both potentials over the momentum range < p < 3 
GeV of relevance for the present calculation. The gaussian has a smaller depth and longer 
range than the OBE potential. Thus the f(p) that determines R embraces a much smaller 
range of kaon momenta for the Barnes' potential than it does for the OBE potential, and 
so gives an R less different from unity in the first case. Moreover, by equating the binding 
energy of the Barnes potential to that of the OBE potential without changing its range one 
finds that it predicts a 1X77) = 0.44 keV, very close to the OBE result. Similarly, the OBE 
potential with its binding set equal to that of the Barnes' potential leads to 1X77) = 0.38 
keV that in turn lies close to the Barnes' result. Thus the calculated two-photon width 
depends mainly on the binding energy of the decaying state. This is a general feature of 
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the weak binding limit k 2 << b 2 . For then it follows from Eqs. (fl9|) . ff2T]) and (fl2l) that 
Tgs ~ V ;2 (0) ~ depends on both the binding energy and potential range parameters 
k and (26) _1 , while the factor R of Eq. ({TBI) varies mainly with the range like R ~ 1/6. 
This is so because its integrand behaves like 1/b 2 times a function of p that suppresses the 
momentum integration beyond p ~ 26. Consequently the range parameter cancels out upon 
forming the product ip 2 (0)R 2 , leading to the behavior 1^(77) ~ k ~ v /(— e). 

The results are summarized in Fig. [2] and Table I respectively. The conclusion is thus 
that there is little to distinguish between these two potential models for describing the 
KK molecule as far the two-photon decay is concerned, once the momentum distribution 
of relative motion of the interacting pair is taken into account. Both lead to very similar 
two-photon decay widths that are in fact in semi-quantitative agreement with experiment. 



5. Comparison with tt-k decay of /o(980) 

The calculated value of r(7T7r) « 50 MeV for the fo — * 2tt decay width already quoted 
in Section 2 for the OBE potential model using the quasi-static approximation falls well 
within the rather wide range of experimental values for r(7T7r) of 40 to 100 MeV reported 
in the literature [t|. On the other hand the gaussian model scales down this result by 
ip 2 (0)/tfj 2 (0) ~ 1/6 to ~ 8.3 MeV that falls way below experiment. It is therefore important 
to check whether non-static effects can also lead to significant changes in either of these 
estimates. 

We assume as before [4J that the transition amplitude for the two-pion decay of the / (980) 
proceeds via K* vector meson t-channel exchange j2l|. Then the photons are replaced by 
pions and the intermediate kaons by K*'s in the first two diagrams of Fig. 1. Note that 
both the direct and crossed diagrams contribute since the the pions are identical bosons in 
the isospin basis. The imaginary part of this combination leads to the analog of Eq. (j4]) 
for pions with the tensor scattering amplitude replaced by a scalar amplitude M n7T . 
Then Eq. ( ITUI) for pions reads Q nw (p) = 87ig 2 nKK 9 n7T (p) where gl nK R is an effective coupling 
constant that determines the free KK — > ttti annihilation amplitude. Its value is not needed 
for the present discussion. Omitting the calculational details, one finds 

Mri = (-4T3-)[ W, ln A/|.. + (p -^- 1 J (16) 

in units Mk, that simplifies to 

1 2p 5 + p z — 2p J 

if we use the estimates Mk* ~ 2Mk and p„ ~ Mk for the K* mass and pion threshold 
momentum in order to compare more directly with 9 11 {p). The function 9 7X . n {p) is shown in 
Fig. 2. Replacing 6* 77 (p) by either form for Q^ip) in Eq. (FiT]) . one obtains a modification 
factor of R^-k ^1.1 for both the OBE and gaussian potentials, an insignificant change. One 
can understand this result without any calculation. It comes about because the momentum 
scale over which 8 nn (p) varies is fixed by the mass of the exchanged boson, in this case 
the K*, that is nearly a factor two more massive than the exchanged kaon that fixes the 
momentum scale for # 77 (p). This means that the value of the integral determining R^ 
becomes almost the same as the normalization integral for f(p) since O^ip) hardly varies 
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at all over the momentum range of f(p). In contrast with the 27 decay problem, the quasi- 
static approximation introduces an insignificant error into the calculation of the 2ir decay 
width. The entries for r(7T7r) in Table I are thus essentially unchanged from their quasi-static 
values. 

6. Discussion and conclusions 

We have shown that calculating the annihilation width of a decaying bound system in 
the quasi-static approximation can be seriously in error when the momentum ranges of 
variation of the annihilation amplitude and the momentum transform of the bound state 
wave function are similar, and have given a revised formulation by including electromagnetic 
corrections in the kernel of the Bethe-Salpeter equation to cover this case. The resulting 
changes in the calculated two-photon annihilation widths for / (980) — > 77 in the molecular 
picture are considerable for one boson exchange model, less so for the Barnes' potential 
model. In fact the predictions for these two competing potentials are brought into semi- 
quantitative agreement with each other and experiment as shown in Table 1. In contrast 
the calculated two-pion annihilation width for /o(980) — > mi is essentially unaffected by 
non-static corrections, but has completely different values for the two models that favor the 
one boson exchange model. 

Since the two-photon decay width in the OBE model has been decreased by nearly an 
order of magnitude by the non-static corrections to already agree semi-quantitatively with 
experiment, the speculation in Ref. jij of important admixtures of, for example, q 2 q 2 states 
in the pure KK ground state in order to reproduce the two-photon width based on the 
quasi-static approximation falls away: The one boson exchange model reproduces the mass, 
as well as reasonable two-pion and two-photon decay widths for / (980) without any further 
assumptions once non-static effects are included in the electromagnetic sector, thereby giving 
additional support to the suggestion that this meson is predominantly a KK molecule. 
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8. Appendix: Equivalent Bargmann Potentials 

For calculational simplicity we replace both the expression for Vobe that appears in 
Eq. (J2j) as well as Barnes' gaussian potential by their Bargmann equivalents that have 
the same respective scattering lengths and effective ranges. These are given by the two- 
parameter potential [13] 




(18) 
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with a = — [1 — \Jl — 2r / ao] / r and 6 = [1 + ^1 — 2r /ao]/r in order to reproduce the 
same scattering length ao and effective range r of the original potential. 

If a = —k < the Bargmann potential in Eq. (fl8l) has a single bound s— state at energy 
e = —k 2 /Mk- The normalized bound state eigenfunction belonging to this eigenvalue is 
given by 



u(r) |-2k(6 2 — a 



r 



2\ 



47T 



1/2 



We note that for weak binding in a short range potential, i.e. ao >> ro, the Bargmann 
potential parameters below Eq. (1181) reduce to a ~ — 1/on and b « 2/ro approximately. 
In this limit one retrieves the standard universal forms [16| e ~ — 1/MkOq and np(r) ~ 
(27rao) _1//2 exp(— r/ao) for the binding energy and radial wavefunction at large distances 
r >> r . 

We also require the Fourier transform 0(p) of '(/'(r) defined by 

m = J -§^Avy p ' r ( 20 ) 

This only depends on the magnitude of the momentum p — |p| since ip(r) in Eq. ffT9|) is 
spherically symmetric. One finds 



/Ajr roc 
d 3 r^(r)e~ ip - r = — ip(0)lm / dr u{r)e ipr 
p Jo 

= y^'^RTTrt 2Fl[1 ' 2 ' 2+ " ; 2 + ^'}' " = 26- ! l 

4, P^oo (21) 

in terms of the hypergeometric function 2-^1 («, A 7; 2) [13]. This result follows immediately 
upon using the substitution x = 1 — exp (— 26r) to transform the radial integral into an 
integral representation of the hypergeometric function. 
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TABLE I: Summary of the calculated values for the one-boson exchange and gaussian potentials. 
The numbers quoted below have been obtained from their Bargmann equivalent potentials as 
described in the text. The experimental values from various sources have also been listed. 

Source ^(0) GeV 3 / 2 ^(77) keV R~ r/ £(77) keV r(7T7r) MeV 

One boson exchange potential [4] 0.091 5.59 0.287 0.46^;^ ~ 50 

Gaussian potential [8j 0.036 0.86 0.632 0.34 ~ 8.3 

Particle Data Group [7J 0.3llg;?f 40 to 100 

Fermilab E 791 Collaboration [18] 44 ± 4 
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FIG. 1: The five diagrams that contribute to order o? to the electromagnetic part of interaction 
kernel in the Bethe-Salpeter equation for K + K~ scattering. Solid lines: charged kaons, broken 
lines: photons. 
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FIG. 2: Variation of the angular integrals 9(p) given by Eqs. (|10p and (|16p with the kaon CM 
momentum for 27 versus 2tt decay of the /o(980) (left hand scale). Also shown are the corresponding 
normalized Fourier transforms f(p) = p 2 (/)(p)/(27r 2 '0(O)) of Eq. f) 1 2 1) for the one boson exchange 
(OBE) and Barnes' gaussian potentials respectively (right hand scale). 
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FIG. 3: Comparison of the numerically generated normalized radial bound state wave functions 
for the one-boson exchange and Barnes' gaussian potentials (solid curves) respectively, with their 
analytic counterparts given by Eq. (|19p for their equivalent Bargmann potentials (broken curves). 
The numerical and analytic wave functions essentially coincide in the case of the Barnes' potential. 
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